Isospin mixing and the continuum coupling in weakly bound nuclei 
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The isospin breaking effects due to the Coulomb interaction in weakly-bound nuclei are studied 
using the Gamow Shell Model, a complex-energy configuration interaction approach which simulta- 
neously takes into account many-body correlations between valence nucleons and continuum effects. 
We investigate the near-threshold behavior of one-nucleon spectroscopic factors and the structure of 
wave functions along an isomultiplet. Illustrative calculations are carried out for the T—l isobaric 
triplet. By using a shell-model Hamiltonian consisting of an isoscalar nuclear interaction and the 
Coulomb term, we demonstrate that for weakly bound or unbound systems the structure of isobaric 
analog states varies within the isotriplet and impacts the energy dependence of spectroscopic factors. 
We discuss the partial dynamical isospin symmetry present in isospin-stretched systems, in spite of 
the Coulomb interaction that gives rise to large mirror symmetry breaking effects. 

PACS numbers: 21.10.Sf, 21.60.Cs, 24.10.Cn, 21.10.Jx 



I. INTRODUCTION 

The charge independence of nuclear force gives rise to 
isospin symmetry [l|, 0] and the formalism of isotopic 
spin has proven to be a very powerful concept in nuclear 
physics [3]. While useful, isospin symmetry is not per- 
fectly conserved. On the hadronic level, isospin is weakly 
violated due to the difference in the masses of the up and 
down quarks 04^] • The main source of isospin breaking 
in atomic nuclei lies, however, in the electromagnetic in- 
teraction 0- 

The members of a nuclear isomultiplet, in particular 
mirror nuclei, provide a unique playground for study- 
ing isospin physics. The invariance under rotations in 
isospin space implies that energies of excited states in 
an isomultiplet should be identical; the deviations are 
usually attributed to the Coulomb force [!, l8j-fll1| . How- 
ever, for nuclear states close to, or above the reaction 
thresholds, the isospin breaking can be modified by the 
coupling to the particle continuum. Here, a spectacular 
example is the Thomas-Ehrman (TE) effect [l2Hl4j that 
occurs when one of the mirror states is unstable against 
particle emission due to a large asymmetry between pro- 
ton and neutron emission thresholds. The resulting TE 
energy shifts strongly depend on the angular momentum 
content of the nuclear state and can be fairly large for 
low partial waves [H, [l6[ . 

The TE effect has also a direct consequence for the 
structure of mirror wave functions [l7H20l | . Indeed, for 
near-threshold states, the configuration mixing involving 
scattering states strongly depends on (i) positions of par- 
ticle emission thresholds in mirror systems (the binding 



energy effect) [2lj], and (ii) different asymptotic behavior 
of neutron and proton wave functions. The latter leads 
to the universal behavior of cross sections [22|, [23| and 
spectroscopic factors (SFs) 0, [H[ in the vicinity of a 
reaction threshold. 

Recently, SFs and asymptotic normalization coeffi- 
cients have been discussed in mirror systems within clus- 
ter approaches (l9l l20j . and strong mirror symmetry- 
breaking in mirror SFs has been predicted. The main fo- 
cus of this work is on the isospin mixing and mirror sym- 
metry breaking in the isobaric analog states (IAS) of light 
nuclei. We show how the different asymptotic behav- 
ior within an isomultiplet and the isospin-nonconserving 
(INC) Coulomb interaction impact wave functions of 
IASs and resulting SFs. Our theoretical framework is 
the complex-energy continuum shell model, the Gamow 
Shell Model (GSM) [H-i|. GSM is a configuration- 
interaction approach with a single-particle (s.p.) ba- 
sis given by the Berggren ensemble [3(| which consists 
of Gamow (bound and resonance) states and the non- 
resonant scattering continuum. 

This paper is organized as follows. Section [TTI presents 
the details of the GSM calculations, with a particular 
focus on the treatment of the Coulomb potential and 
the recoil term. SFs in IASs are discussed in Sec. Mil 
Therein, we study the dependence of SFs on the posi- 
tion of one- and two-particle thresholds. Our calculations 
are performed for prototypical T — l isotriplet consist- 
ing of ,r = 0+ and 2+' IASs in l6 He', 6 Li, and l6 Be\ 
To remove the binding energy effect, we assume identi- 
cal ln/lp emission thresholds. In this way, we isolate 
the effect of the continuum coupling on isospin mixing, 
and study it in the vicinity of proton and neutron drip 
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lines. The results for 6 He, 6 Li, and 6 Be are discussed in 
Sec. IIVI by considering experimental and predicted one- 
particle thresholds. We point out that the conservation 
of isospin in the low-lying states of 6 Be can be explained 
in terms of partial dynamical isospin symmetry present 
in the GSM wave functions of this isospin-aligned system. 
Finally, the conclusions are contained in Sec. [V] 



II. THE MODEL 

The GSM Hamiltonian is diagonalized in the many- 
body Slater determinants spanned upon the Berggren 
s.p. basis. The many-body resonant states of GSM obey 
the generalized variational principle 



31| : they are ob- 



tained using the generalized Davidson procedure that has 
been developed analogously to the generalized Lanczos 
procedure in the context of GSM (see Refs. [H, [29| for 
details). 

We assume in the following that the nucleus can be 
described as a system of n w valence protons or n v va- 
lence neutrons evolving around a closed core. Since our 
discussion concerns the isobaric triplet 6 He- 6 Li- 6 Be, we 
take 4 He as a core. Consequently, the nuclei 5 He and 5 Li 
can be considered as one-particle systems, and 6 He, 6 Li, 
and 6 Be as two-particle systems. In calculations involv- 
ing 5 He and 6 He, the s.p. basis is generated by a Woods- 
Saxon (WS) potential with the radius i?o=2fm, diffusc- 
ness c£=0.65fm, and spin-orbit strength V4o=7.5MeV. 
The depth of the central potential Vb has been varied to 
move the binding energy of a one-neutron system, 5 He 
(i.e., the one-neutron threshold). For Vo=47MeV (the 
" 5 He" parameter set), this potential reproduces energies 
and widths of experimental 3/2j~ and l/2]~ resonances in 
5 He. 

The GSM results should be free from spurious center- 
of-mass (CM) motion. To cope with this problem in our 
GSM approach, we adopt a system of intrinsic nucleon- 
core coordinates inspired by the Cluster Orbital Shell 
Model (COSM) [H|3|. In the COSM coordinates, the 
translationally-invariant GSM Hamiltonian can be writ- 
ten as: 
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where /i is the reduced mass of the nucleon+core system, 
Ui is the one-body WS potential representing the field of 
the core, Vij is the two-body residual interaction between 
valence nucleons, and the two-body term A~ 1 p i pj, with 
A c being the mass of the core, takes into account the 
recoil of the active nucleons. 

The modified finite-range surface Gaussian interaction 
(MSGI) used in this study is a variant of the finite-range 
surface Gaussian interaction (SGI) [Hj]. In order to dis- 
cuss the motivation behind MSGI, we begin with the def- 



inition of the two-body residual interaction SGI: 



V (J,T) exp 



ri - r 2 

Vi 



5(n +r 2 - 2R ) 



1=0 



V (J, T) J2 exp ( — L ^- 1 ) + r 2 ~ 2i? ) 



x f(2£+l)j e [^. )Yt(h)-Y e (r 2 ). 



(2) 



where is the interaction range; Vo( J, T) is the strength 
of the interaction, which depends on the total angular 
momentum J and isospin T; Rq is the radius of the one- 
body Woods-Saxon potential; and f = r/r. 

The contact term represented by the Dirac delta func- 
tion in Eq. ^ generates unwanted divergences in mo- 
mentum space analogous to those present for zero-range 
interactions. To rectify this problem, we replace the ra- 
dial form factors of the multipole expansion of SGI by 
separable terms, chosen independently of £ for simplicity. 
With this choice, the modified interaction MSGI reads: 
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where F(R , 



I + exp 



2Ro - r_p 



(with 



Tf = lfm and fip = 0.05 fm) is a Fermi function which 
makes MSGI practically vanish at r > 2i?o- 

The surface character of MSGI is incorporated through 
the Gaussians centered at Rq. Due to the separability of 
the radial form factors and the presence of the radial 
Fermi cut-off, two-body radial matrix elements of MSGI 
are products of one-dimensional integrals which are non- 
zero only for < r < 2Rq; hence, they are as easy to 
calculate as the radial integrals of SGI [28j]. The range 
of MSGI is fixed at fj,i = 1 fm. The coupling constants 
Vq( J, T) are adjusted to the binding energies ground state 
(g.s.) and first 2 + state of 6 He and 6 Be. It is important to 
point out that the two-body nuclear GSM interaction of 
Eq. ^ is isoscalar by construction. That is, in our work, 
we do not address the question of INC nuclear forces. 

The valence space for neutrons and protons consists 
of all partial waves of angular momentum £ — 0, 1, and 
2. Consequently, the orbital angular momentum cut-off 
in Eq. ([3]) is 4nax=2. The p 3 / 2 wave functions include 
a OP3/2 resonant state and p 3 / 2 non-resonant scattering 
states along a complex contour enclosing the reso- 
nance in the complex fc-plane. For the remaining partial 



waves, i.e. 
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we take the non- 



resonant contour along the real- A; axis (the broad Opi/2 
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resonant state plays a negligible role in the g.s. wave func- 
tion of 6 He and 6 Be). For all contours, the maximal mo- 
mentum value is fc max = 4fm _1 . The contours have been 
discretized with up to 80 points. 

In calculations for systems having valence protons, one 
has to consider explicitly the Coulomb interaction. For 
5 Li, it is represented by a one-body Coulomb potential of 
4 He. In principle, one could approximate it by a Coulomb 
potential of a uniformly charged sphere of radius Rq. 
However, such a potential is inconvenient to use because 
of its non-analytic behavior at Rq. Therefore, we use the 
dilatation-analytic form of the Coulomb potential Uc 
[33 - l36l |. generated by a Gaussian proton density: 

tf<*>(r) = Ze 2 ( 4 ) 
r 

In the above equation, v c — 4Ro/(3y/n), where Rq is 
the radius of the WS potential, and Z is the number of 
protons of the target, e.g., Z=2 for the "proton + 4 Hc 
core" system. The above choice of i?o assures that the 
Coulomb potential given by Eq. Q and the uniformly 
charged sphere potential are equal at r—0. 

The nucleus 6 Be has two valence protons outside the 
4 He core. Consequently, the two-body Coulomb interac- 
tion V c has to be considered. Unfortunately, the calcu- 
lation of two-body matrix elements of V c in a basis gen- 
erated by the one-body part of the GSM Hamiltonian is 
impractical because of difficulties associated with com- 
puting two-dimensional integrals with the complex scal- 
ing method for resonant and scattering basis states. A 
more practical procedure can be developed if one notices 

(2) 

that at large distances the Coulomb term Uc + V c must 
behave as U^(r) ~ 3 e 2 r~ x . Consequently, since uP 
is additive in Z, one can rewrite the Coulomb interaction 
in the 6 Be Hamiltonian as C/ c (3) + (V c - C/ c (1) ). The short- 
range character of the operator V c — Uc suggests using 
the method which consists of expanding two-body oper- 
ators in a truncated basis of harmonic oscillator (HO) 
states [37j : 

v c - up ~ \V C -UW]W 

= P N (V C -UP)P N 

N 

= £ \ap)(ap\V c -UP\-y5)(-y5\ (5) 

where Greek letters label HO states, N is the number 
of HO states used in a given partial wave, and P/v is a 
projector: 

N 

P N = £>/3)<a/3|. (6) 

a/3 

To justify the approximation stated in Eq. ([5]), let us 
consider a normalizable two-body eigenstate l^). |^) can 
be either bound or resonant because resonant states be- 
come integrable when complex scaling is applied to radial 



coordinates [38j. In this case, |\E r ) can be expanded in the 
HO basis used in Eq. ©. According to Eq. @, P/vl*) -> 
|\&) when N — > +oo. Hence, the matrix elements of the 
operator [V c — f/i ]^ involving two-body normalizable 
states converge to those of V c — Uc when N — > +oo. 
i.e. <*/|[^-tf c (1) ] (JV) l*i) -> (9 f \V e -UP\*i). Thelat- 
ter equality is independent of the basis used to expand 
and l^ 7 /). In particular, one can use the Berggren 
basis for this purpose. The short-range character of the 
operator V c - uP implies that (Vf\[V c - [/ C (1) ] (A ° 
should converge rapidly with N. This argument can 
be easily generalized for many-body wave functions with 
more than two particles. 

The matrix elements in Eq. ([5]) can be calculated 
efficiently using the Brody-Moshinsky transformation. 
The computation of one-body overlap integrals between 
Berggren basis and HO states is straightforward, as these 
always converge along the real axis due to the Gaussian 
tail of HO states; hence, no complex scaling is needed. 
The recoil term in Eq. ((!} can be treated in the same 
way as the Coulomb interaction, i.e., by expanding pi in 
a HO basis [37J. The attained precision of calculations 
on energies and widths is better than 0.2 keV for calcula- 
tions without recoil and Coulomb terms, and it is around 
1 keV for the full GSM scheme. 

It has to be noted that because our model involves a 
core, our treatment of the Coulomb interaction is not 
exact. In particular, we neglect the contribution to the 
exchange term arising from the core protons. We also ig- 
nore other known charge- symmetry breaking electromag- 
netic terms such as the Coulomb spin-orbit interaction. 

III. SPECTROSCOPIC FACTORS IN ISOBARIC 
ANALOG STATES 

The SF in the GSM framework is given by the real part 
of the squared norm S 2 of the overlap inte gra l between 
initial and final state in the reaction channel [24, 25|. The 
imaginary part of S* 2 , which is an uncertainty of lZe(S 2 ), 
vanishes if both states in nuclei A and A-\ are bound. 
Using a decomposition of the s.p. channel in the 

complete Berggren basis, one obtains: 

52= |i<*f Hpwa-i 1 ) 3 (?) 

where a^(B) is a creation operator associated with a 
s.p. basis state and the tilde symbol above bra vec- 
tors signifies that the complex conjugation arising in the 
dual space affects only the angular part and leaves the 
radial part unchanged. Since Eq. ([7]) involves summa- 
tion over all discrete Gamow states and integration over 
all scattering states along the complex contour, the final 
result is independent of the s.p. basis assumed. This fea- 
ture is crucial for loosely bound states and near-threshold 
resonances, where the coupling to the non-resonant con- 
tinuum can no longer be neglected. Indeed, the contribu- 
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tion of the scattering continuum to SFs can be as large 
as 25% in such cases [24l [25|. 

In the context of this study, the direct use of Eq. ([7]) is 
impractical when assessing effects related to the configu- 
ration mixing. Indeed, because of the presence of reduced 
matrix elements, S 2 ^= 1 in the absence of many-body 
correlations, and its value depends on j, Ja-i, and J a- 
Hence, we choose to renormalize S 2 by dividing it by 
the extreme single-particle value (obtained by neglecting 
two-body interactions). Within this convention, S 2 = 1 
if configuration mixing is absent. 

The SFs for the two-neutron ( 6 He) and two-proton 
( 6 Be) g.s. configurations considered in our work corre- 
spond to the [ 5 He(g.s.)<g)^p3/2] and [ 5 Li(g.s.)(g>7rp 3 /2] 
channels, respectively. For the T = 1 IASs in 6 Li, 
we consider two channels: [ 5 He(g.s.) <g> np 3 / 2 ]° and 
[ 5 Li(g.s.)®^ 3/2 ] 0+ . 



A. Stability of HO expansion 

The quality of the HO expansion of Eq. (JSJ has been 
numerically checked for both the Coulomb interaction 
and recoil term. Figure [TJ displays the convergence with 
respect to the number of HO states used in the expansion 
for the total energy, width, and SF (real and imaginary 
part) of g.s. configurations in 6 Be and 6 He. For 6 Be, 
the results obtained by assuming the inert core (no re- 
coil, i.e., A c = +oo in Eq. (JTJ) are also presented. It 
is seen that with nine HO states per partial wave, one 
obtains excellent convergence for both energies and wave 
functions, the latter being represented by SFs. For the 
complex energy, the associated numerical error is of the 
order of 2 keV, and this is well below other theoretical 
uncertainties of the model. 



B. Threshold dependence of spectroscopic factors 
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FIG. 1: The numerical check of the HO expansion ([5]) for the 
total energy (top), width (middle), and spectroscopic factor 
(bottom) for g.s. configurations of 6 Be (filled symbols) and 
6 He (open squares). The GSM predictions are plotted as a 
number of HO states used in the expansion. The HO length 
is 6=2 fm [H. 



The GSM SFs for 6 He and 6 Be are shown in the left 
column of Fig. [5J The results are plotted as a function of 
one-nucleon separation energy (Si n for 6 He and Si p for 
6 Be) for three different values of the one-particle thresh- 
old energy Et (i.e., negative of one-nucleon separation 
energy) in one-nucleon systems: 5 He and 5 Li. 

For the bound A=5 systems (Et — —1.5 MeV), the 
SFs in 6 He and 6 Be are different in the whole range of 
separation energies considered. The difference of the SFs 
reach the maximum at the one-nucleon emission thresh- 
old. As the separation energy increases (both nuclei be- 
come more particle-bound) both SFs slowly approach the 
value of one, as expected from simple shell-model consid- 
erations (24|. A characteristic irregularity in the £ = 1 
neutron SF at the neutron emission threshold of 6 He is 
the Wigner cusp. The cusp is absent in the mirror sys- 
tem 6 Be as a result of the different asymptotic behavior 
of the proton wave function [22j ■ 



The energy dependence of SFs changes if the A=5 sys- 
tem happens to be at the particle emission threshold 
(E T =0) or is unbound (E T =0.5 MeV). In both situa- 
tions, a significant difference of SFs in mirror states is 
seen in particle stable (positive Si n or S\ p ) A— 6 systems. 
One may also notice that the Wigner cusp disappears al- 
together if 5 He becomes unbound (cf. Et=0.5 MeV vari- 
ant in Fig. [2]). It is interesting to notice that SFs can be 
greater than 1 if the state of the A — l system is particle 
unstable. This unusual situation (see, e.g., IZe(S) in 6 He 
for Et — +0.5 MeV) is discussed below. 

The imaginary part of the expectation value of an op- 
erator in a resonant state can be interpreted as the uncer- 
tainty in the determination of this expectation value due 
to the possibility of decay during the measuring process 
[39M42| . Figure [5] (right column) shows the uncertainty 
Im(S) of SFs displayed in Fig. [5] (left column). The un- 
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FIG. 2: The spectroscopic factor S = v5 2 , i.e., the real and 
imaginary parts of the square root of the overlap integral 

corresponding to ( 6 He(g.s.)| [ 5 He(g.s.) ® vps,/2\ 0+ ) (solid line) 
and ( 6 Be(g.s.)|[ 5 Li(g.s.)<8>7rp 3/2 ] ) (dashed line) as a function 
of one-nucleon separation energy (Sin for He and Sip for 
6 Be) for three different values of one-particle threshold energy 
Et = —1.5, 0, and 0.5 MeV in A=5 systems (indicated at the 
right of each panel). 




n 


- I I I ; 

\ r \ 


0.04 




0.08 


/ -1.5 J 

;nm|mm|mii|imi: 


0.02 


•\ : 





X \ 


:" h 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 


0.04 







r 0.5 1 


0.04 


, , , i , , , , i , , , , i , , , , : 



-0.1 -0.05 0.05 0.1 -0.1 -0.05 0.05 0.1 

S ln orS lp (UeY) 

FIG. 3: Same as in Fig. [2] except for ( 6 Li(T = l)|[ 5 Li(g.s.) <g> 

^3/ 2 ]° + ) (solid line) and ( 6 Li(T = l)|[ 5 He(g.s.) ® 7rp 3/2 ] 0+ ) 
(dashed line) as a function of one-nucleon separation energy 
(Sin for the former and Si p for the latter) for three different 
values of one-particle threshold energy Et = —1.5, 0, and 
0.5 MeV in A=5 systems. The proton and neutron threshold 
energies are assumed to be identical. 



certainty vanishes if the wave functions in both A=6 and 
A=5 systems are bound with respect to the particle emis- 
sion. Note that in Fig. [2] the appearance of lZe(S)> 1 
cannot be fully explained by the Im(S) plot. Indeed, 
for 6 He at £ r =0.5MeV and 0.05 < S ln < 0.1 MeV, 
TZe(S)> 1 corresponds to \Im(S)\ ~ 0. 

For the J 71 ' = + IAS in 6 Li, we consider two different 
SFs for the (£,j) — P3/2 channel, associated with adding 
a proton to 5 He or a neutron to 5 Li (see Fig. [3]). They 
are plotted in Fig. [3] as a function of one-proton (or one- 
neutron) separation energy. The channel wave functions 
|[ 5 He(g.s.) <8>7rp 3/2 ] 0+ ) and | [ 5 Li(g.s.) <g> vp 3/2 ] 0+ ) are ob- 
viously not orthogonal, as they both share the dominant 
|[ 4 He(g.s.) ® 7T 0p 3 / 2 ® v 0p 3 / 2 } 0+ ) component. The two 
considered SFs for 6 Li differ only by continuum couplings 
induced in the proton and neutron channels. Comparing 
Fig. [2] and Fig. [3J one can see p 3 / 2 proton and neutron 
SFs factors are very similar in both cases. Still, slight 
differences are present. For instance, at Et = 0.5 MeV, 
small irregularities seen in 6 He SFs are absent in the 
neutron SF for 6 Li. A close inspection of proton SFs 
for 6 Li reveals the presence of threshold cusps at zero 
separation, absent on the 6 Be case. This effect can be 
explained in terms of the channel coupling, or flux con- 
servation [24l.l25l|. Indeed, since in our model calculations 
both proton and neutron channels open at threshold en- 
ergy, the coupling between proton and neutron channels 



can generate non-analyticities in proton SFs, even though 
Wigner estimates for proton cross sections are analytical 
at the threshold energy. 

To study the sensitivity of results to the CM treatment, 
we carried out a set of calculations assuming the inert 
core (no recoil). The results are practically identical to 
those of Figs. [2] and [3J The only noticeable difference 
is the absence of a small fluctuation at S , i„«-0.05 MeV 
seen in the real and imaginary parts of SF for 6 He. 

Another reason for the occurrence of lZe(S)> 1 in some 
cases, is the interplay between a final state = 
^"j^-i-R (t ne many-body resonance) and states of the 
non-resonant scattering continuum {^^l^ c } with ener- 
gies close to the resonance energy [39( . The contributions 
to SFs S c in 6 He and 6 Be coming from the non-resonant 
continuum are shown in Figs. 0] and [5] together with con- 
tributions from resonant states. 

In all situations, the contribution of the Gamow res- 
onance to lZe(S 2 ) is dominant. It is interesting to note 
that the impact of the non-resonant continuum does de- 
pend on Et and S\ p / S\ n . For Et < and Si p /S\ n > 0, 
i.e., for A—h and A— 6 bound ground states, the non- 
resonant continuum contribution is basically negligible. 
This is also the case for 6 Be when Et < and Si p /S\ n < 
0, i.e., for a bound 5 Li but an unbound 6 Be. However, 
when either the A=5 g.s. is unbound (Et > 0) or Et < 
and Sip/ Si n < for 6 He, i.e., when 6 He is unbound 
with respect to a bound 5 He, the non-resonant contin- 
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FIG. 4: The distribution of the real part of TLe(Sl) in 6 He 
(solid line) and 8 Be (dashed line) with respect to the 3/2~ 
scattering states of the He and Li systems, ordered ac- 
cording to their real-fc value. The calculations were per- 
formed for three different values of one-particle threshold en- 
ergy S T =-1.5MeV (top), (middle), and 0.5 MeV (bottom) 
in y4=5 systems. The arrows indicate contributions from res- 
onant states. In order to facilitate presentation, the pole con- 
tributions were multiplied by a scaling factor of 0.05. 



uum plays a significant role in both real and imaginary 
parts. In particular, when Et > and S\ p / S\ n > 0, the 
contribution from the non-resonant continuum lZe(S 2 ) 
becomes negative, which translates into a value of lZe(S 2 ) 
that exceeds one. One can also see that Tm(S^) is com- 
parable to lZe(S 2 ), even though this does not occur every 
time lZe(S 2 )> 1. The lesson learned from this discussion 
is that the SF obtained by considering the many-body 
resonance only may often be a poor approximation to 
the total SF, which can contain appreciable non-resonant 
contributions. 



IV. ISOSPIN MIXING IN HE, 6 LI, AND BE 

So far, we have discussed prototypical T—l multiplet 
' 6 He',' 6 Li', and ' 6 Be' with equal proton and neutron sep- 
aration energies to study the effect of different asymp- 
totic behavior on the configuration mixing in the vicin- 
ity of one-nucleon thresholds. In a realistic situation, 
however, particle emission thresholds change within the 
isotriplet due to the Coulomb interaction. To assess this 
effect, we shall now apply the GSM to describe spec- 
tra and SFs for the + ground states and first excited 
2+ states of 6 Hc and 6 Be, and the IASs in 6 Li. In 
calculations involving 5 He and 6 He, we use the " 5 He" 
WS parameter set and the MSGI interaction with the 
strengths: V (J = 0,T = 1) = -15.193 MeV-fm 3 , 



FIG. 5: Same as in Fig. 21 except for the imaginary part of 
Tm{Sl). The arrows indicate contributions from Ops/2 reso- 
nant states if those are not negligible. In order to facilitate 
presentation, the pole contributions were multiplied by the 
scaling factors indicated in the figure. 



V (J = 2,T = 1) = -12.505 MeV-fm 3 . For 6 He, this 
Hamiltonian yields: E a + = -0.974 MeV, E 2 + = +0.823 
MeV, and T 2 + = +89 keV. The experimental values are 



very close: E^* p ^ 



-0.973 MeV, E. 



(cxp) 



and T 



(cxp) 



-n-0.824MeV, 

2+ — -rll3 keV. All binding energies are given 
relative to the binding energy of the 4 He core. 

In the case of 6 Li and 6 Be, we carry out calculations 
in two variants. In variant VI, we take the same WS 
potential as for the He isotopes. Here, isospin is ex- 
plicitly broken by the one-body Coulomb potential and 
the two-body Coulomb interaction between valence pro- 
tons. In variant V2, the depth of the WS potential has 
been changed to 47.563 MeV, in order to obtain an over- 
all agreement for the binding energies and widths of 3/2^ 
and 1/2^ resonances in 5 Li and the + g.s. of 6 Be. The 
readjustment of the one-body potential in V2 is supposed 
to account for the impact of the missing Coulomb terms, 
see discussion at the end of Sec. [Til In both variants, 
MSGI strengths are the same as in the He calculation. 

The predicted g.s. energy of 6 Be, (E + =1.653 MeV, 
r +=41keV) in VI and (E + =1.371 MeV, r 0+ =14keV) 
in V2, is close to experiment: (i?g+ xp ^=1.371 MeV, 

r ( e + x P) =9 2keV). For the first 2+ state, we ob- 
tain: (_E 2+ =2.887 MeV, r 2+ =0.986 MeV) in VI and 
(E 2 + =2.679 MeV, r 2+ =0.804 MeV) in V2. The experi- 
mental energy is (^+ xp) =3.041 MeV, r^ c + xp) = 1.16 MeV). 

Turning to the + IAS of 6 Li, the predicted energy 
is (E + =0.0866 MeV, r 0+ =8.85-10- 3 keV) in VI and 
(E + =-0.0706 MeV, r 0+ =9.13-10- 3 keV) in V2. This 
is fairly close to the experimental value (E^ p ^=- 
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TABLE I: Squared GSM amplitudes of the J" r =0 + IASs of the isotriplet 6 He, 6 Li, and 6 Be. The symbols SI and S2 indicate 
configurations with one and two particles in the non-resonant continuum, respectively. The results for 6 He assuming a rigid 
4 He core are shown in the third column. 



(c fe ) 2 


b He 


6 He (rig. core) 


6 Be (VI) 


6 Be (V2) 


6 Li (VI) 


6 Li (V2) 


(UP3/2) 


U. ( OU— lU.byz 


u. ( yo— iu. i 6i 


1 nnn \(\ o a q 
l.UyU — lU.^4o 


1 1 H7 ;A ICQ 


u.yy4— lu.oo / 


u.y4y— u.oi4 








-0.115+i0.218 


-0.143+i0.255 


-0.084+i0.226 


-0.050+i0.244 


(S'l)^p 3/2 


0.243+i0.619 


0.244+i0.668 






0.066+i0.308 


0.0797+i0.314 


(S2) Sl/2 


0.009+i0.0 


0.0+iO.O 


0.022+i0.0 


0.023+i0.004 


0.011+iO.O 


0.010+iO.O 


(S2) Pl/2 


0.012+i0.0 


0.013+i0.0 


0.008+i0.001 


0.009-i0.0 


0.011+iO.O 


0.012+i0.0 


(S2) P3/2 


-0.049+i0.074 


-0.063+i0.065 


-0.030+i0.029 


-0.028+i0.034 


-0.033+i0.054 


-0.033+i0.055 


(S2) d3/2 


0.002+i0.0 


0.001+iO.O 


0.002+i0.0 


0.002-i0.0 


0.002+i0.0 


0.002+i0.0 


(S2) d5/2 


0.032+i0.0 


0.006+i0.0 


0.025-i0.0 


0.031-i0.04 


0.031+i0.0 


0.031+i0.0 



TABLE II: Same as in Tabled] except for the first 2 + state. 



{Gkf 


6 He 


6 He (rig. core) 


6 Be (VI) 


6 Be (V2) 


6 Li (VI) 


6 Li (V2) 


(Op 3/2 ) 2 




1.132+i0.006 


1.149-i0.022 


0.977-i0.023 


0.987-i0.0267 


1.036-i0.024 


1.049-i0.023 


(Sl)n Pl/2 








-0.004-i0.001 


-0.001+iO.OOl 


0.0-iO.O 


0.0-iO.O 










-0.003+i0.022 


-0.011+i0.027 


-0.001+iO.OOl 


-0.007+i0.0 


(5 , l)y Pl ^ 2 




0.0-iO.OOl 


0.003-i0.002 






0.0-iO.O 


0.0-iO.O 


(5 , l)^ P3/2 




-0.142-i0.009 


-0.147+i0.016 






-0.0492+i0.021 


-0.056+i0.019 


(52)„ Pl/2 


*~P3/2 






0.001-iO.O 


0.001-iO.O 






(S2) vpi ^ 2 


"P3/2 


0.001+iO.O 


0.001+iO.OOl 










( , S'2)p3 /2 




-0.004+i0.006 


-0.005+i0.006 


0.003+i0.006 


0.002+i0.003 


-0.004+i0.008 


-0.004+i0.008 


(S , 2)^d 3/2 


?rd 5/2 






0.001-iO.OOl 


0.001-iO.OOl 






(S'2) J/(i3/2 


" d 5/2 


0.001-iO.O 


0.0+iO.O 










(S2) WI j 3/2 


" d 5/2 










0.001-iO.O 


0.001-iO.O 


(S2) v d 3/2 


7rd 5/2 










0.001-iO.O 


0.001-iO.O 


(S2) d5/2 




0.010-i0.002 


0.0+iO.O 


0.022-i0.006 


0.020-i0.007 


0.015-i0.005 


0.014-i0.003 



0.136 MeV, rj£ xpJ =8.2eV). For the 2+ IAS in 6 Li, we 
obtain (E 2 + =1.667 MeV, T 2 + =0.404 MeV) in VI and 
{E 2 + =1.569 MeV, T 2+ =0.329 MeV) in V2. Both variants 
are in a very reasonable agreement with experimental en- 
ergy: i?2+ xp ^ =1.667 MeV, r£ xp) =0.541 MeV. 

The corresponding g.s. SFs for 6 He and 6 Be (in 
VI) are S" 2 =0.87-i0.383 and 1.015-i0.147, respectively, 
while for 6 Li, they are 1.061— iO. 280 for 7rp 3 / 2 and 
0.911— i0.361 for vp 3j2 . For the 2f state, the SFs 
are S' 2 =1.061+i0.0011 for 6 He, 0.973-i0.0142 for 6 Be, 
0.987-i3.26-10~ 3 for 6 Li (7rp 3/2 ), and 1.034-i0.0235 for 
the 6 Li (vp 3 / 2 ). In spite of the fact that both real 
and imaginary energies in VI and V2 are slightly dif- 
ferent, the SFs for 6 Be in V2 are very close to those 
obtained in VI. Namely, S 2 =1.015-i0.177 for the 
g.s. and S' 2 =0.978-i0.016 for the 2+ state. The V2 val- 
ues of SFs in the 0+ state of 6 Li arc S' 2 =1.028-i0.300 
(7rp 3/2 ) and S 2 =0.898-i0.369 {vp 3/2 ), while for the T=l 
2+ state they are: 5 2 =0.993-i3.190-10" 3 (irp 3/2 ) and 
5 2 =1.043-i0.0224. This can be seen from Table U by 
comparing the corresponding GSM wave function ampli- 



tudes for 6 Be (columns 4 and 5) and 6 Li (columns 6 and 
7). 

While their mean values differ by about 15%, consid- 
ering large imaginary parts, the SFs predicted for the 
jtt = q+ jAgg of the isotriplet, agree within calculated un- 
certainty. However, by examining the GSM wave func- 
tion amplitudes displayed in Table HI one notes that SFs, 
being integrated measures, do not tell the whole story. 
The main effect of the Coulomb interaction is the change 
in distribution of the (0p 3 / 2 ) 2 and (Sl) P3/2 g.s. com- 
ponents, the latter involving one particle in the non- 
resonant p 3 / 2 continuum. As a result, a rather differ- 
ent interference pattern between the resonant 0p 3 / 2 state 
and the non-resonant continuum is predicted for 6 He and 
6 Be, and between resonant and non-resonant states of a 
different type (proton or neutron) in 6 Li. 

For the 2 + IASs, a meaningful comparison of SFs can 
be done as they have small imaginary parts. This is a 
consequence of the smaller configuration mixing induced 
by the nuclear interaction. Indeed, as shown in Table [TT1 
the structure of 2 + states is dominated by the resonant 
(0p 3 / 2 ) 2 component. Here we conclude that GSM pre- 
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diets a mirror symmetry-breaking in SFs of the order of 
5%. 

To assess the impact of the recoil term on our findings, 
we carried out calculations in which the recoil of the core 
is ignored. In this case, the coupling constants refitted to 
the data are V (J = 0,T = 1) = -18.237 MeVfm 3 and 
V (J = 2,T = 1) = -14.942 MeVfm 3 , while the depth 
of the proton WS potential is now 47.5 MeV. Without 
recoil, energy observables are very similar to those ob- 
tained in full calculations. Namely, for 6 He, only the 
width of the first excited state differs by a few keV, as 
it becomes T 2 + — +84 keV. The energy and width of 
the 6 Be g.s. in V2 remain the same as with recoil, while 
there appears a small change for the first excited state of 
6 Be: E 2 + = +2.702 MeV and r 2 + = +0.755 MeV. For 
6 Li, the energy of the + state differs by a few keV in 
VI and around 20 keV in V2, while the width remains 
practically unchanged. For the 2 + state in 6 Li, changes 
are of the order of tens of keV. 

The changes in SFs due to recoil are small as well. 
To show it explicitly, in Table U we compare the GSM 
amplitudes of the ground-state wave function of 6 He in 
the COSM variant (second column) and assuming the 
rigid 4 He core (third column). The main effect of re- 
coil is to slightly redistribute partial wave occupations, 
in particular the {d^/2) 2 contribution. For instance, for 
the 6 He g.s., the sum of the square of amplitudes be- 
longing to the {d^/2) 2 channel is 6T0~ 3 without recoil, 
while it is 3.2- 10 -2 with the full treatment of recoil. For 
Be, not shown in Table |H these numbers in V2 trans- 
late to 4.2-10 -3 and 3.1T0 -2 , respectively. There is also 
a small increase of amplitudes in other continuum chan- 
nels, e.g., (s!/ 2 ) 2 but those wave function components 
are very small. 

Another way of assessing the degree of isospin mixing 
is by inspecting the structure of IAS within the isomul- 
tiplet. To this end, we carried out calculations for the 
isotriplet 6 He, 6 Li, and 6 Be in Vl+COSM using the com- 
mon neutron s.p. basis of 6 He. In this way, the isospin 
operator 

T~ =^aj jT2= _ 1/2 (B)ai 3Tz=1/2 (B) (8) 

is properly defined !43]. The numerical error due to the 
use of neutron s.p. basis on the ground state energy 
of 6 Be is very small: it is about 20 keV for the real 
energy and 5 keV for the with, and this accuracy is more 
than sufficient for the purpose of our IAS analysis. The 
isobaric analogs of the T=l states in 6 He are given by: 

| 6 Li,IAS) = — f-| 6 He), (9a) 

| 6 Be,IAS) = i(T-) 2 | 6 He). (9b) 

The IAS content of a GSM state can be obtained by cal- 
culating its overlap with the state ©. For the + state of 
6 Li, the squared overlap is ( 6 Li| 6 Li, IAS) 2 =0.995. This 



indicates that the lowest + state in 6 Li is indeed an 
excellent isobaric analog of 6 He g.s. Indeed, the corre- 
sponding average isospin value [28J : 



-1 + a/1 + 4(*|T 2 |*) 
T m = V — 2 (10) 

is T Q11 =0.9994. 

For the ground state of 6 Be, ( 6 Be| 6 Be, IAS) 2 =0.951- 
i0.050, i.e., the mean value of the squared amplitude 
exhibits a reduction with respect to the perfect isospin 
invariance. This result is consistent with the large dif- 
ference between GSM wave functions of 6 He and 6 Be: 
a significant component of the 6 Be g.s. wave function 
corresponds to a non-resonant continuum of 6 He. Inter- 
estingly, the total isospin of 6 Be states is perfectly con- 
served in our GSM space. Indeed, having two valence 
protons, wave functions of 6 Be are completely aligned 
in isospace, regardless of the strength of Coulomb inter- 
action. The isospin breaking in 6 Be can only happen 
through core polarization effects, i.e., core-breaking ex- 
citations @, Sj, |45|. Since 4 He is a very rigid core, one 
expects a fairly pure isospin in the low-lying states of 
6 Be. 

A similar situation is expected for any isospin-aligned 
shell-model state corresponding to a semi-magic nucleus 
having Z va \ valence protons (or N va \ neutrons). If one 
disregards core-breaking effects, such a state has pure 
isospin T = Z va \/2 (or T = N va \/2), in spite of the pres- 
ence of INC interactions that manifestly break isospin. 
This is a nice example of a more general phenomenon 
called partial dynamical symmetry, i.e., a symmetry that 
is obeyed by a subset of eigenstates, but is not shared by 
the Hamiltonian [H,II3- We note that while f 2 and f z 
are preserved in the isospin-aligned states, this is not the 
case for T ± operators connecting 6 Be with 6 Li and 6 Li 
with 6 He, that are affected by isospin mixing. 

V. CONCLUSIONS 

There are several sources of isospin (and mirror) 
symmetry-breaking in atomic nuclei. Probably the most 
elusive are consequences of the threshold effe ct 12211 and 
the Coulomb- nuclear interference effect |12hll4| . The 
open quantum system formulation of the GSM makes 
it possible to address the question of the continuum- 
induced isospin symmetry-breaking in a comprehensive 
and non-perturbative way, in terms of the configuration 
mixing involving bound and unbound states. 

As compared to previous GSM studies, present calcula- 
tions are based on a newly developed finite-range residual 
interaction MSGI. The Coulomb interaction and recoil 
term are treated by means of the HO expansion tech- 
nique. The stability of this expansion has been numeri- 
cally checked with a very encouraging result: with only 
nine HO states per partial wave, one obtains excellent 
convergence for both energies and wave functions. 
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To study the sensitivity of results to the CM treatment, 
we carried out two sets of calculations: one in COSM 
coordinates in which the core recoil is treated exactly and 
another one assuming no recoil. We find that the results 
of both variants are very close for both energies and SFs; 
hence, the details of CM treatment do not impact the 
conclusions of our work. 

We have shown that the energy dependence of SFs 
of mirror nuclei is different. Realistic estimates for the 
isotriplet 6 He and 6 Be yield an effect in SFs of the 2 + 
state which is in a range of several percent. This is consis- 
tent with results of recent cluster- model studies [2(| ■ 
For the + configuration, the situation is different. Here, 
the mean values of SFs differ by about 16% and a dif- 
ferent interference pattern between the resonant Qp%/2 
components and the non-resonant p 3 / 2 continuum is pre- 
dicted. However, due to appreciable imaginary parts, 
hence large uncertainty, g.s. SFs in 6 He and 6 Be, and SFs 
for the + analog state in 6 Li, calculated in GSM do not 
offer a clear measure of the mirror symmetry-breaking. 
The behavior of SFs in 6 Li follows that predicted for 6 He 
and 6 Be. Interestingly, proton spectroscopic factors show 
the presence of threshold anomalies due to the strong 
coupling with the neutron channel. 

Due to the partial dynamical isospin symmetry present 
in the GSM wave functions of 6 Be, the low-lying states in 
this isospin-stretched (T=l, T z =—1) system are expected 
to show very weak isospin breaking effects. This is in 



spite of Coulomb interaction present in this nucleus. For 
the T z =0 member of the isotriplet, 6 Li, the isospin sym- 
metry is explicitly broken in the GSM space as a result 
of mixing between T=0 and T=l states but the result- 
ing mixing is very weak. We thus conclude that the large 
mirror symmetry breaking effects seen in binding energies 
and SFs of the isotriplet are related to T ± components 
rather than the total isospin. 

In summary, the coupling to the non-resonant con- 
tinuum can give rise to isospin and mirror symmetry- 
breaking effects that are configuration dependent. Ex- 
planations of mirror symmetry breaking based on the 
traditional close quantum system formulation of the nu- 
clear shell model sometimes invoke INC nuclear effective 
interactions [HI EH . We would like to point out that any 
attempt to extract such interactions from spectroscopic 
data should first account for the coupling to the many- 
body continuum in the presence of isospin-conserving nu- 
clear forces. If neglected, or not treated carefully, the 
continuum effects can alter results of such analyses. 
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